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A new operator for pion photo- and electroproduction has been developed for nuclear applica- 
tions at photon equivalent energies up to 1 GeV. The model contains Born terms, vector mesons 
and nucleon resonances (P33(1232), Pii(1440), L>i3(1520), 5'ii(1535), Fi5(1680), and L»33(1700)). 
The resonance contributions are included taking into account unitarity to provide the correct phases 
of the pion photoproduction multipoles. The dependence of electromagnetic resonance vertices 
QQ I is described with appropriate form factors in the electromagnetic helicity amplitudes. Within this 

. model we have obtained good agreement with the experimental data for pion photo- and electropro- 

^\ ' duction on the nucleon for both differential cross sections and polarization observables. The model 

can be used as a starting point to predict and analyze forthcoming data. 
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Pion photo- and electroproduction is presently one of the main sources of our information on the structure of 
nucleons and nuclei. With the advent of the new generation of high intensity, high duty-factor electron accelerators 
\ as Jefferson Lab (Newport News), MAMI (Mainz) and ELSA (Bonn), as well as modern laser backscattering facilities 
' as LEGS (Brookhaven) and GRAAL (Grenoble), this field reaches a new level of promise. 

' The elementary amplitude of pion photo- and electroproduction on free nucleons is one of the main ingredients of 
C3 . the analysis of these reactions for nuclei. It has been the subject of extensive theoretical and experimental studies, 
and over the past 30 years a series of models was developed for photon energies from threshold up to 500 MeV (see for 
example Refs. and references therein). Attempts to extend these models to energies up to around 1 GeV, by use 
of effective Lagrangians and coupled channels approaches, have been presented in Refs. |^,|^ and [||. Unfortunately, 
'^T' most of the recent models are too sophisticated and complicated and difficult to handle for nuclear applications. 

A simple and effective model for nuclear calculations was developed in the late 70's by Blomqvist and Laget It 
^ provided an adequate description of then available experimental data up to the first resonance region. However, in 
P5 . the mean time modern electron accelerators have provided a host of new high precision data. Beams of high current 
^ ' and high duty- factor have reduced the statistical errors to the order of a few percent, and the new data have provided 
• ^ us with high quality multipole analyses of pion photoproduction from VPI and Mainz jl^] . 

With such advances it is certainly appropriate to develop new models for the description of pion photo- and electro- 
H [ production. With respect to nuclear applications, these models should have a simple and physically transparent form. 
- - ' On the other hand they should satisfy gauge invariance and unitarity, and reproduce the existing experimental data, 
not only in the first but also in the second and third resonance regions, for both pion photo- and electroproduction, 
which are the subject of the current interest. It is therefore the aim of this work to develop a model having these 
properties. 

As a starting point wc will use the prescriptions of the isobar model [ p2|p^ ] , assuming that resonance contributions in 
the relevant multipoles have Breit-Wigner forms. The dependence of the jNN* vertices will be determined via the 
corresponding helicity amplitudes or quark multipole moments [ p^ . The nonresonant contributions will be described 
using standard Born terms with a mixed pesudovector-pseudosclar ttNN coupling and vector meson exchange. The 
final amplitude will be unitarized by extending the procedure developed by Olsson ||l^ and Laget jl^ to the case of 
virtual photons and higher resonances. Finally we demonstrate that such an extremely economical model provides 
a good description for individual multipoles as well as differential cross sections and polarization observables. We 
believe that the developed model can be used not only for nuclear applications, but also as a starting point to predict 
and analyze forthcoming data for pion photo- and electroproduction on proton and neutron targets. 



'Permanent address: Laboratory of Theoretical Physics, JINK Dubna, Head Post Office Box 79, SU-IOIOOO Moscow, Russia. 
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II. FORMALISM 



In the following we will briefly summarize the most important expressions which clarify the conventions used in our 
paper. For details we refer to, e.g. Ref. 

In accordance with Ref. 0| we decompose the current operator in the ttN cm frame into the standard CGLN 
amplitudes F, [llira 



rriN 



ia-Fi + {a ■ q){cr X k)F2 + iq{(T ■ k)F3 



+iq{(T ■ q)F4 + ik{cr ■ fe)F5 + ik{(T ■ q)Fe 



(1) 



47rW 

rriN 



i{(7 ■ q)Fr + i{(T ■ k)Fs 



k J 



where totv and cr are the mass and spin operator of the nucleon respectively, W is the total energy of the ttN system, 
k = k/ \ k \ and q = q/ \ q \ are the unit vectors for the photon and pion momenta respectively, and a = a — {a ■ k)k 
represents a vector with purely transverse components. The eight amplitudes Fi, -Fs a-re complex functions of three 
independent variables, e.g. the total energy the pion angle 07r, and the four-momentum squared of the virtual 
photon, = —up' > 0. Due to current conservation we obtain the two relations | fc | F5 = wFg and \ k\FQ = uoFi. 
Therefore only six amplitudes are independent. In the present work we use i^i, -Fg- 

For the analysis of the experimental data and also in order to study individual baryon resonances, pion photo- and 
electroproduction amplitudes are usually expressed in terms of three types of multipolcs, electric {Ei±), magnetic 
(Af/±), and longitudinal {Li±) ones, with pion angular momentum I and total angular momentum j = I ± 1/2. They 
are defined by a multipole decomposition of the amplitudes Fi, 



Fi = J2^ilMi+ + Ei+)Pl+, + [{1 + l)Af,_ + Ei^]P;_,}, 

l>0 

F2 = Y.l^l + l)Mi+ + lMi_]P!, 
i>i 

Fi = Y,{Mi+-Ei+-Mi^-Ei^)Pl\ (2) 

l>2 

where P/ are the derivatives of the Legendre polynomials. Note that in the literature the longitudinal transitions are 
often described by Si± multipoles, which correspond to the multipole decomposition of the amplitudes P7 and Eg,. 
They are connected with the longitudinal ones by Si± =| k \ Li±/uj. 

From total isospin conservation in the pion-nucleon system it follows that the amplitudes Fi (or the multipoles) can 
be expressed in terms of three independent isospin amplitudes [|l^. These are A^^^ for the isoscalar photon, and for 
the isovector photon the two amplitudes A^^/^^ and A^^/^^ for the TriV system with total isospin 7=1/2 and / = 3/2 
respectively. However, it is also useful to define the proton pA^^/"^^ and neutron „A(^/^) amplitudes with total isospin 
1/2, 



^^(1/2) ^ ^(0) ^ 1^(1/2) ^(1/2) ^ ^(0) _ 1^(1/2) 

3 ' 3 



(3) 



With this convention the physical amplitudes for the four physical pion photo- and electroproduction processes are 



A(7*p^ n7r+) = V2 
A(7*n pTr~) = \pi 



4(1/2) _ 1^(3/2) 

3 
3 



A[^*n ^ nTfO) = -„A(V2) ^ 2^(3/2) 

3 



(4) 
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All the observables for these processes may be expressed in terms of the amplitudes Fi. In general there are 
16 different polarization observables for the reactions with real photons. In pion electroproduction we have four 
additional observables due to the longitudinal amplitudes F5 and Fq and 16 observables due to longitudinal-transverse 
interference, giving a total number of 36 polarization observables. In the literature there occur many different 
definitions for these observables. In our work we shall mainly follow the convention of Ref. ||l^. In the case of 
comparison with experimental data presented in different conventions, we will give special comments. 



III. BACKGROUND 



A. Born and vector mesons exchange terms 



As pointed out in the introduction, the main goal of our present work is to develop a simple model for nuclear 
applications, which should be consistent with the recent data for pion photo- and electroproduction on nucleons at 
photon energies (or equivalent energies) up to 1 GeV. This energy range covers the first and second, and touches the 
third resonance regions. The base line for a correct description of the resonance contributions is, of course, a reliable 
description of the nonresonant part of the amplitude (nonrcsonant background). Traditionally this part is described 
by evaluation of the Feynman diagrams derived from an effective Lagrangian density. For the electromagnetic 'yNN 
and 77r7r vertices the structure is well defined. 



jO^nn = -dip 



7^^^Fr"(g") -I- ^9M''F|'"(Q2) 



V', (5) 



C-,^^ ^ e[{d^T^)^ X T,]^A^ F^{Q^) , (6) 

where is the electromagnetic vector potential, ^ and tt are the nucleon and pion field operators, respectively. In 
Eqs. (5-6) we have included explicitly dependent proton (Ff 2)' neutron (^"2) ^^d pion (F^) electromagnetic 
form factors. In the case of real photons the form factors are normalized to ff (0) = ^^(0) = 1, i^i"(0) = 0, ^I'(O) = 
Hp — 1.79 and F2{Q) = Kn — —1.97. For virtual photons we express the nucleon form factors in terms of the Sachs 
form factors by the standard dipole form, and assume that F^((3^) — Ff{Q^) — Fi{Q'^) which is the simplest way to 
preserve gauge invariance. In the same way we treat the axial form factor Fa{Q^)- 

The description of the hadronic ttNN system is a more sophisticated part of the theory of pion photo- and electro- 
production. In this case there are two possibilities for the construction of the interaction Lagrangian, the pseudoscalar 
(PS), 

^Inn = igi'l5T ■ 1/'^ , (7) 

and the pseudovector (PV), 

CllN^-—^l5l^'r-d''n^, (8) 

ttNN coupling, where g'^/'in = 14.28 and //m^r — g/2mM- At low pion energies the PV coupling is to be preferred, 
because it fulfills PCAC and is consistent with low energy theorems (LET) and chiral perturbation theory to leading 
order. However, the PV model cannot be renormalized in the usual way, and this produces a problem at high energies. 
On the other hand, the renormalizable PS model leads to a better description at the higher photon energies. In our 
work we consider both of these schemes by a gradual transition between them. 

The final expressions for the CGLN amplitudes Fi, ...,^6, obtained using the effective Lagrangians (5-8) are well- 
known (see, for example Ref. fl^ ). We will refer to them as Born term contributions. They are the dominant part 
of our background. The other part is related to vector meson exchange contributions. In general they are much 
smaller, but as we will show below they are quite important for some multipoles. The effective Lagrangians for uj and 
p exchange are Q 

£ = e^e^,p„{d^'A'')^,dP{5aLo'' + p",)Fv{Q^), (9) 

CvNN = [gviit. + 1^^/^-^") (^^ + ^ • P^)V' , (10) 

where oj and p are the uj and p meson field operators respectively, and Ay is the radiative coupling determined by 
y — > 777 decay. The dependence of C-^t^v is defined using a dipole form factor. 
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In general the values for the strong coupling constants gvi (vector coupling) and gv2 (tensor coupling) in Eq. (10), 
are not well determined. In various analyses they vary in the ranges ||^^po| of 8 ^ g^i ^ 20, 1.8 ^ gpi ^ 3.2, — 1 ^ 
51^2/5^1 ^ and 4.3 ^ .9^2 /Spi ^ 6.6. In the present work we take them as free parameters to be varied within these 
ranges. The off-shell behaviour of the vertex functions is described by hadronic monopole form factors. 



B. Nonresonant multipoles 

From the considerations above we find that the following ingredients for the construction of the background are 
not well defined: 1) the type of the ttNN coupling and 2) the coupling constants for the vector meson exchange 
contributions. The best way to fix them is to analyse nonresonant s- and p-wave multipoles. In our energy region 
these are E^^/'^\ m[^J'^\ p,„Mi+^^\ and p.,,E[^/'^\ At photon lab energies E^ < 500 MeV, the multipole p,„£^^?/^^ 
can also be considered as nonresonant. 

First, let us fix the coupling constants of the vector mesons. For this purpose we consider the multipoles which are 

independent of the type of the ttNN coupling. These are p, ,1-^14/^^ and p,n£'i+^\ of which the first one (Afi+) is 
especially sensitive to uj exchange contributions, as illustrated in Fig. 1. The final results for the coupling constants 
and cut-off parameters Ay are given in Table 1. 

The real parts of the nonresonant multipoles Eq^ and A/i_ in the isospin 3/2 channel are more appropriate to fix 
the optimum parametrization for the ttNN coupling. In Fig. 2 we have depicted their energy dependence obtained 
with PS (dotted curves) and PV (dashed curves) couplings. In the threshold region the best description of recent 
analyses from VPI [0 and Mainz can be obtained within a PV model. As pointed out before, this coupling has 
also to be preferred on general grounds, because it reproduces the leading terms of ChPT near threshold. However, at 
larger photon energies {E^ > 600 MeV) the results of the VPI group are in between our calculations with PV and PS 

(3/2) (3/2) 

couplings. In such a situation, in order to describe Eq_^ and Afj^_ multipole in a wide energy range up to E^ = 1 
GeV, the most economical way is to construct a hybrid model (HM) with a mixed type of ttNN coupling. In order 
to give a reasonable threshold behavior, it has to start with pure pseudovector coupling at the lowest energy and will 
develop into pseudoscalar coupling at the highest energies. The simplest effective Lagrangian for such a model may 
be written in the form 

A 2 2 

'—kNN — A 2 :„2'-ttNN + A 2 :„2'-TrNN i \^^) 

where Qq is the asymptotic pion momentum in the nN cm frame which depends only on W and is not an operator 

(3/2) (3/2) 

acting on the pion field. From the analysis of the ' and i?o+ multipoles we have found that the most appropriate 

value for the mixing parameter is A,„ = 450 MeV. We note that for the pion pole terms the pion couples with on-shell 
nucleons only, and HM, PV and PS models are all equivalent. Moreover, only the multipoles Eg^, Mi_, Lo+, and 
Li_ are affected by changing the coupling schemes. 

One of the important peculiarities of the i?o+ multipoles is the relatively large imaginary part, even at low pion 
energies. The origin of this feature is well-known: pion-nucleon rescattering (or final ttN interaction). In order to take 
account of this effect, we shall use a prescription in accordance with unitarity (Fermi- Watson theorem) and K-matrix 
approach, i.e 

E^„2 = E^^liBorn + u,p){l + i tij^) , (13) 

where — [777 exp(i(5^^) — l]/2i is the pion-nuclcon elastic scattering amplitude with the phase shift S^j^ and the 
inelasticity parameter 777 (both taken from the VPI analysis). 

In Fig. 2 (solid curve) we see that the unitarity condition (13) combined with our Lagrangian (12) provides an 

(3/21 

excellent description of the real and imaginary parts of the Eqj^ multipole in a wide energy range. In Fig. 2 we 
illustrate similar results for the case of isospin 1/2 and E^ < 550 MeV. In the next section we will describe the analysis 
of the latter channel for the higher energies in connection with 5*11(1535) resonance excitation. 

Finally we note that in the charged pion channels the multipole i?o+ is not sensitive to different choices of the 
ttNN coupling and contributions of pion rescattering. However, these effects are extremely important in the neutral 
channels for both the Eq^ and Lo+ multipoles. In the latter case we predict a large imaginary part that certainly 
has to be taken into account in future analyses of the (e, e'lr^) reaction, e.g., in the transverse-longitudinal (TL) cross 
section, where the nonresonant /m{io+} multipole interferes with the large resonant Im{Mi^} multipole. 
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IV. RESONANCE CONTRIBUTIONS - REAL PHOTONS 



The background contributions being fixed, we can develop a reliable scheme to study baryon resonances by analyzing 
the relevant multipolcs. However, it is well-known (see e.g. Ref. that even in this case the procedure for the 

extraction of the "bare" resonance contributions is not trivial due to the interference with the background. In 
the present work we generally consider so-called "dressed" resonances which include "bare" resonances and vertex 
corrections due to the interference with the background. We believe that for nuclear applications this is a more 
appropriate way to facilitate the investigation of medium effects. On the other hand, the K-matrix approach asserts 
that at resonance position the contribution from the interference term is small or even vanishing in the case of an 
ideal resonance. Therefore, we expect to get reliable informations about the conventional resonance parameters at 
the position of the resonance. 

First we consider pion photoproduction. For the relevant multipolcs Ai± we describe the resonance contributions 
assuming a Breit-Wigner energy dependence of the form 



Ai±{w) = A± MW) fMw) c^N , (14) 



where /7rjv(W^) is the usual Breit-Wigner factor describing the decay of the N* resonance with total width Ttot, partial 
ttA^- width Tttn and spin j, 



UNiW) 



1 kw rriff T^ 



N 



(2j -I- 1)^ I 9 I w rl 



1/2 



The factor Cj^n is y/3/2 and -l/\/3 for the isospin 3/2 and isospin 1/2 multipolcs respectively, as defined by Eq. 
(3). 

In accordance with Refs. [pp3] the energy dependence of the partial width Tj^n is given by 









' Wr 








W 



r.. = P.TR^^) [j^)^, (16) 

where X is a damping parameter, assumed to be X = 500 McV for all resonances. Tji and are the total width and 
the pion cm momentum at the resonance peak (W — Wr) respectively, and /3,r is the single-pion branching ratio. 

The total width Ttot in Eqs. (14-15) is the sum of F^^v and the "inelastic" width Tin- For the latter one we assume 
dominance of the two-pion decay channels and parametrize the corresponding energy dependence as in Ref. |2^ , 

Ttot = r.AT + F,„ , F„ = (1 - f3^) Tr (^) ^ ( + f \ ^ , (17) 

\qoJ \X^ + q^^J 

where q2Tr is the momentum of the compound (27r) system with mass 2m7r and qo — q2-^ at VF = Wr. Concerning the 
definition of Fi„ , it takes into account the correct energy behavior of the phase space near the three-body threshold. 
We make an exception for the 5*11(1535) resonance, where we also introduce a r\l^ width, similar to Eq. (16) but with 
the mass = 547 MeV. 

The main parameters in the 'yNN* vertex are the electromagnetic amplitudes Ai±, introduced in Ref. |2^. They 
are linear combinations of the usual electromagnetic helicity amplitudes A1/2 and (see, e.g., Eq. (1) in Ref. p2t). 
We parametrize the W dependence of the jNN* vertex beyond the resonance peak with a form factor 



where the damping parameter X is the same as in Eqs. (16-17) and kR = kw at IF = Wr. In order to preserve the 
correct " pseudothreshold" behaviour, mainly given by the Born terms, we introduce a parameter n ^ Ij, with l-y the 
orbital angular momentum of the photon. 

One important ingredient of the Breit-Wigner parametrization (14) is the unitary phase </>. The main role of this 
parameter is to adjust the phase ip of the total multipole (background plus resonance) to the corresponding pion- 
nucleon scattering phase St^n (in accordance with the Fermi- Watson theorem, when the influence of the inelastic 
channels is small) or to the experimentally observed phase. The latter procedure will become necessary when the 
photon energy increases above 500 MeV, the approximate limit of the Fermi- Watson theorem. 
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For the most important A(1232) resonance (m}^^^-* and Ef'_^'^^ multipoles), the influence of the inelastic channels 
is negligibly small up to « 800 MeV. In this region ■0 = <5^Ar with very good accuracy. At the higher energies, 
where the A contribution becomes small, we will use the ansatz fE3| 



ij^liW) = arctan 



'^-riiiW) cos25l{W) 



ril{W) sin2(5f(VF) 



(19) 



The pion-nucleon scattering phases (5/ and inelasticity parameters r]l are taken from the analysis of the VPI group 



(SAID program). Then, applying the unitarization procedure developed by Olsson |15| and Laget |l6|, we can 
determine the phase (l){W). The results of our fit to the Ml^l ' and EY_l ' multipoles are shown on Fig. 3, and the 
corresponding phases 4>{W) are given in the appendix. 

We apply the same procedure to the Pii(1440) resonance and the corresponding multipoles pM'\^J_'^'^ and „Mji^^''. 
However, in this case the influences of the inelastic channels on the pion photoproduction phase become important 
already at E^ > 450 MeV, and they are different in the proton and neutron channels. In this situation we define 
4>{W) by using the results of the VPI analysis. 

In the second and third resonance regions the most important resonances are the 5*11(1535) contributing to the 

Eq]^'^^ multipole, the I?i3(1520) with contributions to the E^i"^^ and m'^I'^^ multipoles, and the J^i5(1680) with 

contributions to the E'^i'^^ and multipoles. In the isospin 3/2 channel the most important resonance is 

the 1)33(1700) which contributes to the i?2_ and multipoles. Unfortunately, for all these multipoles little 

information is available about the phases beyond the resonance peak. Therefore, in these cases we consider a free 
parameter that we take as a constant extracted from the observed ratio between the imaginary and real parts at the 
resonance peak. 

The final results for the more important multipoles and the corresponding values for the resonance parameters are 
presented in Figs. 4-6 and Tables 2-3, respectively. We note that the -fNN* vertex at resonance is described by the 
two helicity amplitudes ^1/2 and ^3/2. These can be easily extracted from the imaginary part of Eq. (14) using the 
values for the electromagnetic amplitudes Ai±, the phases ipR = (j){W = Wr), and Eq. (1) of Ref. For example, 
we find in the case of the A resonance 

Ai/2 - COS^E + M^^'f^ cos0m)/2, 

A3/2 = V3{£fJ^^ cos - A^f/'^ cos 0Af )/2 . (20) 

Our values for the helicity amplitudes are summarized in Table 4. In general they are in the range of the recent 
Particle Data Group (PDG) analysis ||2^], the only exception being the case of the neutron Pn amplitude. 

In Figs. 7-8 we compare our results (solid curves) for differential cross sections, photon (E) and target (T) asym- 
metries with the recent dispersion relation analysis of Mainz ||ll| (dashed curves), and with old and new data from 
Bonn and Mainz in the first resonance region. One can see that our model describes the observables in this region 
at the same level of accuracy as dispersion relations. Moreover, in this energy region our model can be extremely 
simplified by keeping only contributions from ^33(1232) and Pii(1440) resonances, and the calculations with pure PV 
ttNN coupling and with our hybrid model differ very little. 

In Figs. 9-12 we present our results (solid curves) for the differential cross sections and single polarization observables 
at E-y > 450 MeV, where the contributions from the second and third resonance regions become increasingly important. 
Our complete model generally agrees well with the experimental data in the TT~^n channel over a wide energy region 
up to 1 GeV. However, there is one exception, the differential cross section at 700 MeV. From Fig. 4 we can see that 
this is the region where the 5*11(1535) resonance has a maximum contribution due to a cusp effect. In this region 
calculations without the 5ii resonance (dashed curves) describe the data better. In the ir'^p channel we could only 
obtain a satisfactory agreement with the data up to E^ < 900 MeV. Unfortunately, the data sets of different groups 
are not consistent in that channel, and the error bars are much bigger than in the case of the data at lower energies. 

In Fig. 13 we demonstrate the evolution of the energy distribution in the tt+ti channel at backward angles by 
subsequent addition of resonance contributions to the Born terms. First, we see that the large contribution from (non 
unitarized!) Born terms is strongly suppressed by the A resonance. The second interesting feature is related to the 
excitations of the 5ii(1535) and Z)i3(1520) resonances. At 6^^ ~ 120° their contributions are comparable and it is 
difficult to separate them. However, at 9-^ = 180" the contribution of the Z?i3(1520) resonance is very small, because 
the E2~ multipole vanishes at this kinematics. This fact provides a possibility to study the contribution from the 
5*11(1535) resonance and the related cusp effect due to the opening of the channel, clearly separated from 1)13(1520) 
effects. 
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V. RESONANCE CONTRIBUTIONS - VIRTUAL PHOTONS 



The major problem in the extension of om' model to the case of virtual photons is associated with the determination 
F the dependen 
this dependence by 



of the dependence of the amplitudes Ai±{Q'^) at the resonance peak. For the main Al^'^^'' amplitude we define 



)(g2) ^^W2)(o)LMe-7QV^(Q2)^ (21) 

where Fd is the usual dipole form factor. The parameter 7 will be determined later. Note that in Eq. (21) an 
additional dependence appears due to the kinematical factor | k \ /k\Y, with the virtual photon three-momentum 

At present time there is only little experimental information available on the dependence of the small E^j^^'^ 
and S'^'^^^'' multipoles. Following Ref. we therefore assume that the corresponding electromagnetic amplitudes 
have the same dependence as 

^ kw 

S^ff\Q') = £^^f\Q') = ^-^C^^f\Q^) . (23) 

Of course, this assumption will have to be checked whenever more information about the £'i_|_/Mi+ and S'i+/Afi+ 
ratios will be available from the new experiments. 

The extension of the unitarization procedure to the case of pion electroproduction is straightforward for the A- 

(3/2) (3/2) (3/2) 

resonance, because the Fermi- Watson theorem requires that the phases of the total M^_^[ ' , El_^ ' and S{_^ multipoles 
should not depend on Q^. Therefore, the position of the resonance peak does not depend on either. In our model 
this important requirement is satisfied by an appropriate choice of the phase (f) in Eq. (14). We note that after 
the unitarization procedure this phase will depend not only on W but also on Q^. In this way we always obtain 
resonance multipoles with real parts vanishing at the resonance position W — 1232 MeV and with the total phases 
of the multipoles equal to 90". In Fig. 14 we give our predictions for the and multipoles at Q^^O, 0.2, 

and 1.0 (GeV/c)^ which satisfy this requirement. Note that the i?^'^^^'' multipole has a second zero immediately above 
the resonance peak. Its position is dependent, and in accordance with unitarity the real and imaginary parts 
simultaneously change sign. The obtained unitary phases as a function of and W are given in the appendix. 

In Fig. 15 we show our unitarized M^^'^'^ multipole as a function of at W — 1232 MeV. The experimental data 
were extracted from the magnetic G^^ form factor using the relation 



Mir\Q')-'-^J^f^GlAQ'), (24) 
mjv y STexpqR 



where a = 1/137 and Texp — ll5MeV. The best fit of the experimental data was provided with a parameter 7 = 0.24 
(c/GeV)^ in Eq. (21). In Fig. 15 we also show results obtained without additional unitarization for finite Q^, i.e. 
<j) — <t>R =const, independent of Q^. In this case the real part of the M^j^^^'' multipole does not vanish at the resonance 
as soon as is finite (dotted curve). 

Important information about Z?-state components of the nucleon and A wave functions may be extracted from 
the study of the Rem = Ei^/Mi+ and Rsm = 'S'i+/Mi+ ratios in the isospin 3/2 channel. In Fig. 15 we give 
our predictions for the dependence of these ratios at resonance. The ratios are practically constant in the non- 
unitarized case of = 4>ii=C0Tist (dashed curves), because we assumed that the dependence of the £'1+ and 
multipoles is the same as in the case of Mi+. However, this behaviour changes after unitarization, especially for 
Q2 < 0.5 (GeV/c)2). 

Only little information is available on the multipole pM^J^"^^ related to the excitation of the Roper resonance 
Pii(1440). In our analysis as well as in many other ones, we neglect the contributions from longitudinal resonance 
excitations, because of a lack of any reliable data. Moreover, an analysis of the existing exclusive and inclusive 
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electroproduction data |Q showed no indication of an excitation of the Roper resonance at high Q^. This imphes 
a rapid fall-off of the corresponding proton hehcity amplitude A1/2 with . However, such a result would be at 
variance with numerous quark model calculations. 

Following Ref. ji^ we take the dependence of the electromagnetic amplitude pTM^^^'* at resonance {W = 1440 
MeV) in the form 

Ml'\Q') =P -P(-£) > (25) 

where a is the harmonic oscillator parameter of quark model calculations whose value varies from a = 0.229 GeV, 
which reproduces the proton charge radius, to a = 0.410 GeV, which fits the A3/2 helicity amplitude for the Z?i3 
resonance Beyond the resonance peak we introduced an additional dependence in the form factor /^jv using 
the usual prescription, replacing kw ^| ^ | and kn ^ kn =\ k \w=Wri- 

Due to the importance of the inelastic channels, we can no longer apply the Fermi- Watson theorem to the Pn 
resonance as well as to resonances in the second and third resonance regions. As in the case of real photons we 
will assume that the real part of the multipoles vanishes at resonance, i.e. 4'{Wr) — 'n jl independent of . As a 
consequence of Eq. (14), the phase becomes a dependent function. The procedure of finding (l){Q^) at resonance 
we will again call unitarization procedure. 

As may be seen seen in Fig. 16, the application of the unitarization procedure to the pM^]/^'' multipole strongly 
modifies the original dependence defined by Eq. (25), leading to a rapid fall-off with Q^. In the non- unitary 

approach, where = 0ij=const, the dependence is smoother and the real part of the pM^]/^'' multipole has a large 
negative value at resonance. In Fig. 16 we illustrate the effect of unitarization also for the proton helicity amplitude 
Ai/2 which is related to the amplitude pTW^^^^ by 

A/2(Q') ^pM['I^Hq') COS ^{Q^). (26) 

For the harmonic oscillator parameter we have taken the value a = 0.370 GeV. With such a value our unitarization 
procedure works up to Q^=l (GeV/c)^, while smaller values of a limit applicability to smaller momentum Q^. 

The dependence of the pion electroproduction amplitudes in the second and third resonance regions has been 
defined in accordance with Refs. |l|,||,^, where 5ii(1535), Di3(1520), £'33(1700), and ^15(1680) resonances are 
considered the dominant states of the [70, l"]i and [56, 2+]2 super multiplets in an SU{6) (g) 0(3) symmetry scheme. 
Within this model the transitions between the nucleon and its resonances are expressed in terms of the quark electric 
and magnetic multipoles ej^^ and m^"^ , where A is the photon helicity, L and J = L or L -I- 1 are the orbital and total 
angular momenta, depending on the multipolarity of the transitions. As in most other analyses of the experimental 
data, we neglect the contributions from the longitudinal (or Coulomb) transitions assuming that the corresponding 
quark multipoles Cg ^ and ^ vanish. 

Our choice to express the dependence in terms of the quark multipole moments is motivated, first, by their 
simple connection with the amplitudes Ai± listed in Table 5. Second, as we will see below, the quark multipole 
moments exhibit a very simple dependence. The third advantage is that by this way we can obtain both proton 
and neutron amplitudes. As may be seen in Table 5 we have assumed a certain deviation from the predictions of 
the SU (6) symmetry in the case of the neutron amplitudes in order to match with our values for the corresponding 
helicity amplitudes at the photon point, leading to a 10% increase for the £'13(1520), a 10% decrease for the £15(1680) 
and 30% decrease for the S'ii(1535) resonances. 

In terms of the quark multipoles the experimental data have been compiled by Burkert et al. in the form 

el'^el'FDiQlvp), ml' ^rhl' Fd{QIvf) (27) 

for the transition to the [70, l^]i super multiplet describing the excitation of the S'ii(1535), Di3(1520) and £33(1700) 
resonances, and 

ef = gf Qevf Fd [QIvf) , ™f = rhj-' Qevf Fd (QIvf) (28) 

for the transition to the [56,2+]2 super multiplet describing the excitation of the £15(1680) resonance. In order to 
minimize relativistic effects, the dependence is evaluated in the equal velocity frame (EVF) with 



Am^Wn ^m^Wn ' ^^^^^^^^ " (1 + QW/0.71)2 " (29) 
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The dependence of the so-cahed reduced quark multipole moments ef"^ and mf"' is parametrized similarly as in 
Ref. [Q with a shght modification in order to match with the amphtudes Ai± of Table 3 (or the helicity amplitudes 
of Table 4) at the photon point. For the transition to the [70, l^]i super multiplet we obtain 

e}i = 2.56, = im-imQ\yp, 

-0.42 + 5.82 Qlyp, Qlyp ^ 1.0 GeV^ 

5.84 - 0.44 Qlyp, Qlyp > 1.0 GeV' 



(30) 



For the transition to the [56, 2+]i super multiplet the parametrization is 

ef=4.35, mf = 3.823 - 4.172 Q|y^, 
.23 „ / 0.42 + 5.1 g|^^, Qlyp ^ 1.0 GeV2 



"1 \ 6.31 - 0.79 Qlyp, Qlyp > 1.0 GeV^ . ^^^^ 

Note that in Eqs. (30-31) Qevf, e^"' and rh^"^ are given in GeV, e^'^ and rh?"^ are dimensionless. In Figs. 17 and 
18 we demonstrate the consistency of our parametrization with the reduced quark multipole moments and helicity 
amplitudes compiled by Burkert [^,|4j. Note that for the second 5*11(1650) resonance we followed the assumption of 
Ref. H, 

^Sii(i650)(q2) ^ Af)f'\Q') tan 30° , (32) 

where the 30° mixing angle matches our value for the helicity amplitude at the photon point. 

The unitarization procedure for the excitation of the Z3i3(1520), D33(1700) and Fi5(1680) resonances is similar to 
the procedure for the Pii(1440) resonance, i.e. we assume that the total phase tp — 7t/2 atW — Wr is independent 
of . In Fig. 19 we compare the dependence for the more important proton electric and magnetic multipoles 
at resonance to the results of the analyses from DESY, NINA and Bonn, compiled by Foster and Hughes . In 
contrast to the case of the Pn resonance, our unitarization procedure for these multipoles does not seriously modify 
the original dependence of the quark multipole moments, as can be seen by comparing the solid and dashed curves 
in Fig. 19. 

In the following Figs. 20-23 we present some examples for observables of 7r+ and 7r° electroproduction. In general 
the results of our model are in good agreement with the existing data. In particular we note that our results obtained 
without Pii(1440) resonance contribution get closer to the data with increasing Q^. This indicates a strong decrease 
of ^1/2(^11) as function of . We hope that future experimental studies will clarify these problems occurring for the 
Roper resonance contribution. Furthermore a precise L/T separation would be very useful in order to obtain more 
information about the dependence of the Si- multipoles. 

Finally, we comment on the u-channel resonances which are often present in effective Lagrangian approaches. Obvi- 
ously, the main subject of experimental and theoretical studies are the s-channel resonances, while the contributions 
from the tails of the u- channel resonances are just slowly varying backgrounds which are essentially real and affect 
all the multipoles. We assume that the w-channel contributions are effectively included in the unitarization procedure 
and the nonresonant background of our model. More sophisticated methods beyond the isobar model of our work will 
demand an explicit inclusion of the M-channel resonances in a crossing symmetric way. 



VI. CONCLUSION 



We have developed a new operator for pion photo- and electroproduction for applications to reactions on nuclei at 
photon equivalent energies up to 1 GeV. The model contains Born terms, vector mesons and nucleon resonances up 
to the third resonance region (P33(1232), Pii(1440), i:»i3(1520), 511(1535), Fi5(1680), and 1533(1700)). 

For the Born terms we propose an energy dependent superposition of pseudovector and pseudoscalar ttNN couplings. 

This procedure describes the correct energy dependence of the nonresonant multipoles for photon lah energies up to 

3/2 3/2 

— 1 GeV, in particular in the case of the and M^'_ multipoles, and provides an excellent agreement with 

the results of the Mainz multipole analysis at photon energies below 500 MeV and with the VPI analysis at higher 
energies. 

The resonance contributions in the s-channel are parametrized by the standard Breit-Wigner form for the relevant 
multipoles. Unitarity is fulfilled by multiplying the resonance terms with appropriate phases. We find that the 

(3/2) (3/2) 

unitarization procedure developed for pion electroproduction is very important for extracting the Mj^_^ , and 
LfJ_'^^ resonance multipoles, especially for small Q^. In the case of the Pii(1440) resonance, we assume that the phase 
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of the total pAfji"' multipole is independent of Q"^ at resonance which leads to a suppression of the Roper resonance 
with increasing values of Q^. 

The Q"^ dependence of the electromagnetic form factors for the Sn, D13, 1)33, and F15 resonances is expressed in 
terms of quark multipole moments, attributing the first three states to the [70, 1~] and the last state to the [50, 2+] 
multiplets. The assumption of a independent phase at resonance leads only to a small modification of the 
dependence of the hclicity amplitudes in these cases. 

Within our model we obtain good agreement with the existing experimental data for pion photo- and electropro- 
duction on the nucleon. Due to its inherent simplicity, the model is well adopted for predictions and analysis of future 
results on pion electroproduction on nucleons and nuclei. 
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APPENDIX: 

In this appendix, we give the W and dependence for the phases <^ in the Breit-Wigner paramctrization obtained 
after the unitarization procedure. First, let us consider the W dependence at the photon point (5^=0. In this case 
our results for the phases of the magnetic, electric and longitudinal components of the ^33(1232) resonance can be 
parametrized by 

(/)m(0) = 22.130a; - 3.769a;2 + 0.184a;3 , 

0^(0) = 83.336a; - 28.457a;2 + 3.356a;'^ - 0.122x4 , 

,^,^(0) = 43.668.x - 2.872a;2 - l.giOx^ + 0.230a;'* , (Al) 

where x = {W — m„ — mjv)/(100 MeV). All phases are given in degrees. 

For the phases of the proton and neutron components of the magnetic excitation of the Pii(1440) resonance we 
have 

(l)p{0) = -4.661a; + 0.349a;^ - 0.240a;^ , 

(l)n{0) = -22.935a; + 19.935a;2 - 4A56x^ + OMSx^ + 0.039a;^ . (A2) 

The longitudinal excitation of the Roper resonance is neglected in our model. For the other resonances we assumed 

the phases (j) to be constants as listed in Table 3. 
The dependence of (p is parametrized as 

«^) = ^*<°'- <*^' 

The coefficients A and B for the P33(1232) resonance depend on W and are given in Table 6. With good accuracy 
the phases of all other resonances can be taken as constants, identical for proton and neutron. The corresponding 
values for these phases are given in Table 7. Note that for the and resonances we take A = B = 0. 
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mv[MeV] 


\v 


gvi 


gv2/gvi 


Av[GeV] 




782.6 


0.314 


21 


-0.57 


1.2 


p 


769.0 


0.103 


2. 


6.5 


1.5 



TABLE I. Masses and coupling constants for the vector mesons. 



N* 



Wfl[MeV] 



rfl[MeV] 



13^ 



Multipoles 

p(3/2) »^(3/2) 

^(1/2) „(l/2) 



^33(1232) 
Pi 1 (1440) 
I>i3(1520) 
S'ii(1535) 
S'ii(1650) 
Fi5(1680) 

L'33(1700) 



1235 
1440 
1520 
1520 
1690 
1680 
1740 



130 

350 

130 

80 

100 

135 

450 



1.0 

0.70 

0.60 

0.40*' 

0.85 

0.70 

0.15 



p(l/2) 

p(i/2) 

p(l/2) 
p,n-'-'3_ ; p,n-'^^3 
^(3_/2) ^(3_/2) 



(1/2) 



TABLE II. Parameters of nucleon resonances. *' The branching ratios for the eta and two-pion decay channels are 50% and 
10%, respectively. 



AT* 


Ai± 


4>R 


Ai± (proton) 


(neutron) 


n 


^33(1232) 


£.(3/2) 
''1+ 


26 
73 


323 
-17 


323 

-17 


2 
1 


Pii(1440) 




-25 


-78 


66 


1 


L'i3(1520) 




23 

35 


-145 
-68 


149 
23 


2 
4 


S'ii(1535) 
S'ii(1650) 


?(l/2) 

?(l/2) 
P,n<-'0+ 






-67 
-39 


55 
32 


3 
4 


Fi5(1680) 


<?(l/2) 


15 
15 


-64 
-37 


7.7 
22 


3 
4 


1)33(1700) 


£(3/2) 


61 
61 


-240 
38 


-240 
38 


4 
4 



TABLE III. Unitary phases (fiR at the resonance position (in deg) and electromagnetic amplitudes Ai± (in 10 ^ GeV ^^^), 
n is the parameter in Eq. (18). 
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N* 


1 A 7 \ 

^1/2 (proton) 


A 7 VI 

^3/2(protonJ 


1 A 7 ^ 

A.-1 /2 (neutron j 


A 7 \ 

-^3/2 (neutronj 


^33(1232) 


our 


-Loo 


-ZOD 


-loo 


-ZOD 




PDG 


1 Ar\4- ^ 


-ZOo^ D 


1 40-1- 


-ZOo^ D 


Pi 1 (1440) 


our 


- / i 




fin 






PDG 










L'i3(1520) 


our 


-i / 


1 fi/l 




-ioo 




PDG 


9/14- Q 


iDD^ 




1 QQ-I- 1 1 


S'ii(1535) 


our 


D / 




-00 






L LJ\J 


70-1- 1 9 








511(1650) 


our 


39 





-32 







PDG 


53ib 16 




-15ib 21 




Fi5(1680) 


our 


-10 


138 


35 


-41 




PDG 


-15± 6 


133± 12 


29± 10 


-33± 9 


£>:',:', (1700) 


our 


86 


85 


86 


85 




I'DG 


10 15 


85^ 22 


10 15 


S.")^ 22 



TABLE IV. Helicity amplitudes (in W-^GeV'^^^). 



N* 


multiplet 


Ai±Cr cos (j)R 


proton 


neutron 


£'13(1520) 


[70,1-] 


p,nC'2- C R COS (pR 

p,nM z^^^Cr COS 4>R 


-V2e" - m" 


1.1 (^26" + im") 


Sii(1535) 


[70, 1-] 


p,n£'()+ ^^CkCOS 4>r 


-e" + V2m'' 


0.7(e--y|m") 


Pis (1680) 


[56,2+] 


p,n£^i^J^^CRCOS (t>R 
p,nM3^I^^CR COS (j)R 




0.9y^|m- 
0.9 ^fm^^ 


£'33(1700) 


[70, 1-] 


p,«f2- '^'C'flCOS 4>R 
p,nM^l^^CR COS (j)R 


-V2e" + im" 


-V2e" + |m" 



TABLE V. Electromagnetic amplitudes Ai±Cr cos (f>R, where Gr= |Y^3mAr(W^ — m\,). Fourth and fifth column: these 
amplitudes expressed in terms of the quark multipole moments, for proton and neutron respectively. 



multipolos 








W{MeV) 


A 


B 


A 


B 


A 


B 


1100 


1.500 


3.859 


5.377 


27.29 


1.802 


19.64 


1150 


1.227 


2.794 


4.070 


15.65 


1.825 


10.75 


1200 


0.9147 


2.175 


3.036 


10.25 


1.422 


6.739 


1250 


0.8444 


1.762 


2.800 


7.886 


1.393 


4.587 


1300 


0.8320 


1.456 


2.490 


5.252 


1.267 


2.991 


1350 


0.7810 


1.225 


2.046 


3.262 


0.9363 


1.715 


1400 


0.7277 


1.045 


1.775 


2.239 


0.5052 


0.7424 


1450 


0.6815 


0.9035 


1.633 


1.720 


0.3361 


0.3500 


1500 


0.6471 


0.7963 


1.604 


1.500 


0.2979 


0.2100 



TABLE VL The coefficients A and B for the ^33(1232) resonance (in GeV~^). 
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A and B 












A 


9.977 


15.17 


1.563 


-0.070 


-1.024 


B 


2.384 


21.98 


5.202 


3.766 


3.938 



TABLE VII. The coefficients A and B for the Pii(1440), L»i3(1520) and Fi5(1680) resonances (in GeV~'^). 



14 




450 600 750 
E (MeV) 



150 300 450 600 750 900 1050 
E (MeV) 



n 1+ 



~l — ' — ' — T" 




ii 

,*A/\A m- 



I I I I I I I I I I I I I L 



-1 

150 300 450 



F (1/2) 
n 1+ 



'AATA A / 



1- 




600 
E (MeV) 



750 900 1050 



FIG. 1. Nonresonant multipoles p,nM[^^^^ and p,n£'i+ "^^ calculated without (dashed curves) and with (solid curves) vector 
meson contributions. The tvNN coupling constant is taken as (?^/47r — 14.28 and the coupling constants for the vector mesons 
are given in Table 1. The open and full circles are the real and imaginary parts from the Mainz dispersion analysis jl^. The 
full and open triangles are real and imaginary parts from the VPI analysis |l^ . 



150 300 450 



600 
E (MeV) 



750 900 1050 



r{l/2) 




150 300 450 600 750 900 1050 150 
E (MeV) 

5 ' 



250 350 450 

E (MeV) 



550 



a -5 




I 

o 



- > -10 - 




-10 -15 

150 300 450 600 750 900 1050 150 250 350 450 550 

E^(MeV) E^(MeV) 

FIG. 2. Nonresonant Eo+ and Mi_ multipoles calculated with pure pseudovector (dashed curves) and pure pseudoscalar 

(dotted curves) nNN couplings. The solid curves are the results for the real and imaginary parts obtained using the Lagrangian 

(12) and the unitarization ansatz (13). Data points as in Fig. 1. 
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FIG. 4. The multipoles p,„M^i/^' and j,,„B^+^^' calculated with the Pii(1440), 511(1535) and S'ii(1650) resonances. The 

dashed curves for the p,nM^_l'^^ multipoles are the Born + u>,p contributions with our Lagrangian (12). The real (dash-dotted 

curve) and imaginary (dotted curve) parts of the p,ni?o+^^ multipoles are obtained without the S\i resonances using the 
Lagrangian (12) and Eq. (13). Data points as in Fig. 1. 
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, 5. The multipoles E2- and M2- calculated with the _Di3(1520) and 1)33(1700) resonance contributions. The dashed 
are the Born + uj,p contributions. The full and open triangles are the real and imaginary parts from the VPI analysis 
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FIG. 7. Differential cross sections, photon (E) and target (T) asymmetries for p(7,7r"'")n in the first resonance region. The 
solid and dashed curves are the results of our calculations and the Mainz dispersion analysis respectively. Experimental 
data from Refs. 125 
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